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Abstract 

We extend the formalism of the thermodynamic two-time Green's functions to nonextensive 
quantum statistical mechanics. Working in the optimal Lagrangian multipliers representation, the 
g-spectral properties and the methods for a direct calculation of the two-time g-Green's functions 
and the related g-spectral density (q measures the nonextensivity degree) for two generic operators 
are presented in strict analogy with the extensive {q = 1) counterpart. Some emphasis is devoted to 
the nonextensive version of the less known spectral density method whose effectiveness in exploring 
equilibrium and transport properties of a wide variety of systems has been well established in 
conventional classical and quantum many-body physics. To check how both the equations of 
motion and the spectral density methods work to study the g-induced nonextensivity effects in 
nontrivial many-body problems, we focus on the equilibrium properties of a second-quantized 
model for a high-density Bose gas with strong attraction between particles for which exact results 
exist in extensive conditions. Remarkably, the contributions to several thermodynamic quantities 
of the g-induced nonextensivity close to the extensive regime are explicitly calculated in the low- 
temperature regime by overcoming the calculation of the q grand-partition function. 
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I. INTRODUCTION 



The method of thermodynamic Green's functions (GFs) ,Q is a 

powerful tool in ordinary statistical mechanics for exploring the equilibrium and transport 
properties of a large variety of many-body systems. These functions are related to important 
physical quantities and hence their calculation constitutes one of the basic problems of 
extensive thermostatistics. 

Remarkably, the extraordinary effectiveness of the GF technique in quantum many-body 
physics [l, 2, s], li, 0, [g, Q] has stimulated a lot of research activity to extend this successful 



method also to study extensive classical statistical mechanics js, 9,ll0|. In the context of the 



two-time GFs 
and Sadovnikov 



the foundations of the classical formalism were introduced by Bogoliubov 
four decades ago and further systematic developments in this direction 



were performed only many years later j^, 9]. Additionally, the classical counterpart of the 
quantum Matsubara GF framework \^ was achieved in Ref. |lo|. So also the two-time GF 
formalism in extensive classical statistical mechanics can be now considered well established. 
Although not currently used in the literature, it has been successfully employed to study 
;he thermodynamics and the transport properties of several classical many-body systems 
2I, Is, ^, 10] also involving phase transitions and critical phenomena j^. 

Various (also numerical) methods have been developed for the calculation of the two- 
time GFs both in quantum and classical statistical physics . In addition 
to the most known "equations-of-motion method" (EMM), the related "spectral density 
method " (SDM), originally formulated by Kalashnikov and Fradkin [3| within the quantum 
statistical mechanics context, appears to be a very promising nonperturbative approach to 
perform reliable studies of the macroscopic properties of classical and quantum many-particle 
systems avoiding an explicit calculation of the partition function. 

At the present time, the situation does not appear so well established in nonextensive 
statistical mechanics which has attracted an increasing interest since the seminal proposal by 



Tsallis made almost 20 years ago 
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121]. This framework can be regarded as a generalization 



of the Boltzmann-Gibbs statistical mechanics to properly describe macroscopic memory 
effects (e.g., non-Markovian stochastic processes) and, generally speaking, systems exhibiting 

mm 



nonergodic microscopic dynamics 
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Despite its elegant formalism, the Tsallis thermostatistics is affected by some intrinsic 
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difficulties in performing explicit analytical calculations for realistic many-body systems. 
Nevertheless, by resorting to different techniques and approximations as generalizations of 
the extensive ones, the Tsallis theory has been successfully applied to a wide variety of 
systems for which nonextensivity effects are not negligible 13|, [ij, llSj and must be taken 
into account for a proper comparison with experiments. Known theoretical tools have been 
employed and adapted in the Tsallis framework such as, for instance, linear response theory 



16l | , p erturbation and variational methods [17|], path integral 
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Monte Carlo 



2J] and molecular dynamics [2l| techniques and many others. Less attention. 



however, has been devoted to the extension, to the nonextensive many-body world, of the 



well-established thermodynamic GF technique llOj l in ordinary 

quantum and classical statistical mechanics. 

n 

Some years ago the GF method in the Kadanoff-Baym framework [4| was generalized to 
the Tsallis quantum statistical mechanics adopting the second quantized representation for 
many-particle systems 25|, 1261]. In these works, the q GFs for a nonextensive many-body 
system were formally expressed in terms of parametric integrals over the corresponding 
extensive (g = 1) quantities, q denoting the so-called Tsallis parameter which measures the 
nonextensivity degree. Of course, this method may be really useful when the many-body 
problem for the extensive counterpart has been solved but, unfortunately, this is possible 
only in a limited number of situations. In any case, the crucial step to calculate nontrivial 
contour integrals in the complex space constitutes an additional formidable problem which 
would require, generally, further approximations. Thus, many-body methods which allow 
one to perform direct calculations of the ^properties, overcoming the a priori knowledge 
of the related extensive ones, are desirable. Motivated by the conviction that a direct q 
GF method may provide new and effective calculation techniques to deal with nontrivial 



nonextensivity problems, we have recently extended [27|, l28| the Bogoliubov-Sadovnikov 



91] to nonextensive classical statistical 



two-time GF framework j^, and the related SDM 
mechanics working conveniently within the so-called optimal Lagrangian multipliers (OLM) 
representation suggested in Ref. This choice avoids some intrinsic difficulties involved 
in other ones and allows simplified analytical and numerical calculations. In any case, our 
suggestion is quite general and can be extended to different contexts preserving the physical 
content, consistently with the equivalence of the current four versions of the Tsallis statistics 
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The aim of the present article is to extend the two-time GF formahsm to the quantum 
nonextensive statistical mechanics within the OLM representation. Particular emphasis will 
be devoted to the spectral density (SD) and its spectral decomposition due to their relevance 
for practical calculations. 

The q EMM and the q SDM are here presented as powerful tools for a direct calculation of 
ihe two-time q GFs and q SD, respectively. Besides, in parallel with the extensive counterpart 
qI, we outline the key ideas to explore the dispersion relation and the damping of elementary 
and collective excitations in nonextensive many-body systems. Finally, with the aim to show 
how the two-time GF method works when the Tsallis q distribution is involved, we present 
analytical calculations for a nontrivial high-density Bose gas with str ong attraction between 



the particles for which exact results exist in the extensive case |31l . |32| |. The effects of the 
g-induced nonextensivity will be consistently explored in the low-temperature regime using 
both the q EMM and q SDM. 

The article is organized as follows. In Sec. [Tll we present a summary of the basic 
ingredients of the nonextensive quantum statistical mechanics which will be useful for the 
next developments. In particular, we focus on the quantum OLM representation which 
appears more convenient for our purposes. Section IIIII is devoted to the formulation of 
the two-time GF method in the Tsallis quantum statistics, with emphasis on the spectral 
properties. In Sec. llVt we present a nonextensive version of the EMM and SDM. Here, the q 
SDM is properly implemented to offer the possibility to study systematically the dispersion 
relation and the damping of excitations within a unified formalism. Section |V] deals with 
the application of these methods to a model which describes a high density Bose gas with 
strong attraction between the particles. Concluding remarks are drawn in Sec. IVI[ Two 



appendixes close the article. In Appendix A the recently proposed [27|, |28| two-time GF 
method in Tsallis classical statistical mechanics is shortly reviewed with the aim to point 
out the substantial differences between the nonextensive quantum and classical frameworks. 
Appendix B summarizes some mathematical details. 
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II. BASIC INGREDIENTS OF NONEXTENSIVE QUANTUM STATISTICAL 
MECHANICS 



The Tsallis nonextensive thermostatistics constitutes today a new paradigm in the field 
of statistical mechanics. The key problem of the Tsallis framework is to find the appropriate 
von Neumann density operator p which maximizes the generalized q entropy (in units = 
h=l) 

s, - (1) 

subject to appropriate constraints related to the evaluation of the mean or expectation value 
of the observables within the nonextensive scenario. Here Tr{. . .) stands for the usual trace 
operator. 

In the literature, four possible choices have been considered for the evaluation of q ex- 
pectation values {■■■)q: 

I — I 

(i) The original Tsallis proposal [11] 

{A), = TripA), (2) 

where the Hermitian operator A corresponds to a generic observable A. For a system 
with Hamiltonian H, this implies the canonical representation 

p=[(l-q)(a + PH)/q]^, (3) 

in terms of the two Lagrange multipliers a and f3. Unfortunately, this choice involves 
some troubles related to the Lagrange multipher a. 



(ii) The Curado-Tsallis (CT) choice 33] 



{A), = Tr{p^A), (4) 

which yields the canonical result 

P = Z-i[l-(l-g)/5i/]^, (5) 

with 

Z, = Tr[l~{l-q)(3H]^^ . (6) 

This avoids the explicit presence of the multiplier a but has the disadvantage to exhibit 
unnormalized mean values ((l)q 7^ 1). 



(iii) The Tsallis-Mendes-Plastino option [1 



(7) 



implying the canonical solution 

p = z;' 



with 



Zq — Tr 



1-9 



Trip'l 



{H-Uq 



1-9 



(8) 



(9) 



where Ug = {H)q is the q internal energy. Here we have (1)^ = 1 normalized q mean 
values, but troubles occur again in obtaining p due to the presence of Tr (p'^) = {Zq 
in Eqs.([S]) and (self-referential problem). 



(iv) The OLM improvement 29|] which preserves the gentropy ([T]), but replaces the Tsallis- 
Mendes-Plastino-like constraints by "centered" mean values. Essentially, the general 
OLM procedure consists in maximizing the Tsallis generalized entropy ([I]) subject to 
the constraints 

TrP = 1 (10) 

and 

Tr[p'^iA,-{A,)q)] = (J = 1...,M), (11) 

where the generalized mean values {Aj)q = T^ (p'^Aj) /Tr{p'^) for M relevant Aj {j = 
1,...,M) observables are assumed to be known a priori and hence are regarded as 
constraints in the variational approach to the nonextensive thermostatistics. 

In the canonical representation this yields 



p = Z^'[l-{l-q)(3{H-Uq)]^^ 
where the generalized partition function Zq is now given by 

Zq = Tr[l-il-q)PiH-Uq)]^^. 



(12) 



(13) 



As we see, the OLM framework avoids all the inconveniences occurring in the previous 
choices. Besides, in Eqs. f lT2|) and f|T3l) . the Lagrange multiplier /3 does not depend on 
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the partition function and it is identified as the inverse physical temperature [29|, |3C 
[jS = 1/T), consistently with the zeroth law of thermodynamics. 

The extension of the canonical OLM prescription to the grand-canonical ensemble can 
be simply obtained from Eqs. f|T2l) and f|T3l) by replacing H and Uq with Ti = H — fiN 
and Uq = {'H)q = Uq — ^{N)q, where fi is the chemical potential and denotes 
the operator describing the total number of particles in a system. Other ensemble 
representations can be obtained similarly. 

It must be emphasized that the four versions (i)-(iv) of the Tsallis thermostatistics are 
equivalent in the sense that the probability distribution for each of them can be easily 



derived from any other of them by using appropriate transformation rules jSOj. In any case, 
for practical calculations one must select the most convenient version case by case. 

As already mentioned, in the present work we find it convenient to adopt the OLM 
version to formulate the two-time GF method in nonextensive statistical mechanics, but the 
framework can be easily extended to different versions by using the prescriptions given in 
Ref. (30i. 

To close this short review of the basic elements of the Tsallis thermostatistics, we mention 
two relevant features of the nonextensive scenario. 

The first one, tacitly assumed before, is that the basic OLM canonical operator fq = 
l — {l — q)f3{H~Uq) must be positive definite. This means that, for any specific nonextensive 
problem, one must take into account only microstates \n), with H \n) = En\n), satisfying 
the cutoff condition [12] 

< (T^^ , if ? < 1> (14a) 

E^>—^-— + Uq ,ifg>l, (14b) 
(1 -g)/5 

where the quantity — q)f3 + Uq could have, in principle, positive or negative sign. 

The second crucial question concerns the physical origin of the nonextensivity parameter 
q. Unfortunately, little has been definitevely proposed in terms of first principles to explain 
the appearance of the Tsallis thermostatistics in many real situations and to understand the 
physical meaning of the elusive parameter q. We cite below some general ideas which, in 
our opinion, may have a seminal role in future research activity on the subject. 

After several atte.p. fl, a well «„Hed scenario emerged almost seven years ago, 



now called the "finite-heat-bath picture". First Almeida 



34| derived the Tsallis power-law 



probability distribution from first principles assuming exactly constant the heat capacity 
Cb of the thermal bath in contact with the system of interest. Here, the physical meaning 
of q is simply expressed in terms of C^- Specifically, when it is finite, q ^ 1 and one has the 
Tsallis distribution, while for an infinite heat capacity of the bath, q = 1 and one recovers 



the conventional exponential distribution. This picture has been further elaborated 35|, |36| , 
but a substantial progress was performed in Ref. jsG] where, via a model-free derivation of 
the Tsallis statistics (i.e., without resort to the microscopic details of a system and of its 
surrounding), it is shown that the finiteness of the heat capacity of the environment is not 
necessarily due to only the finiteness of its degrees of freedom, as reductively assumed in 
previous treatments. 

A more general point of view was reported in Ref. [37|] where the parameter q was inter- 
preted as a measure of fluctuations of the parameters (as the temperature) in the exponential 



distribution and expressed in terms of the variance of their inverse. The superstatistics 
subsequently developed as a generalization of this idea, consists in a superposition of different 
statistics relevant for driven nonequilibrium systems with complex dynamics in stationary 
states with large fluctuations of intensive parameters (e.g., the temperature, chemical po- 
tential, energy dissipation, etc.) on long time scales. It explains the emergence of the power- 
law statistics for real systems as a result of fluctuations in their environments, supporting 
the Tsallis framework as a special case. In this "fluctuational picture," the nonextensivity 
parameter q is found to be directly related to proper stochastic processes or constraints 
imposed to the systems and, in the absence of fluctuations, the extensive case is consistently 
reproduced. Remarkably, the superstatistics seems to offer a general formalism for treating 
nonequilibrium stationary states of complex systems which exhibit dynamics that can be 
decomposed into several dynamics on different time scales. Some effort has been performed 
recently jsol to legitimate this theory as a true statistical mechanics framework, but several 
aspects remain to be further elaborated and clarified. 

Other, rather fragmentary, proposals to justify the Tsallis statistics exist, but we mention 
here only some of them that, in our opinion, may be relevant in nonextensive many-body 
physics. 

There is a broad consensus of opinion that the nonextensivity may occur also in cases 
where long-range interactions play a relevant role, the gravitational and Coulomb forces 



being important examples with laboratory and astrophysics implications [15]. However, 
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no definitive statistical foundation for long-range interacting systems, such as the self- 
gravitating ones, has still been well established and, in particular, the true physical nature of 
q in terms of the range of the forces and the space dimensions has not been well understood 



yet 



40|. 



Also confi ning traps for interacting bosonic or fermionic systems may induce nonexten- 
sivity effects |4l| due to the interplay between the particle interactions and the trapping 
potential. Moreover, interesting findings about the physical origin and the measurement 
of the parameter q have been obtained recently for complex magnetic systems such as the 
manganite LaQ,jSrQ,^MnO^/ MgO and the amorphous alloy CwgoCoio (42[. Such compounds 
seem to embody three basic ingredients which may induce nonextensivity: long-range in- 
teractions, clusters with fractal shapes, and intrinsic inhomogeneity. With these features in 
mind it has been shown [42] , via theory (at a mean -field approximation level) and scanning 
tunneling spectroscopy measurements, that the magnetic properties of these materials can 
be properly described using the Tsallis thermostatistics, q measuring the competition of the 
intrinsic inhomogeneity and dynamics. 

Finally, for our demonstrative application to a high-density Bose gas performed in a 
next section, we have in mind some of the nonextensivity mechanisms outlined before. For 
instance, we think about a system of bosons in a fluctuating environment ,^ in confining traps, 
or confined to a self-gravitational field as in non-relativistic boson stars [43]. 

We conclude this section noting that it is frequently convenient to use the so-called q 
exponential function 

e^^[l + (l-g)x]^, (15) 

which simplifies sensibly the formalism of the nonextensive statistical mechanics. One can 
verify that this function (a) yields e\ = for g ^ 1, (b) for g > 1 vanishes as a power 
law when x — oo, and (c) diverges at x = —1/ (1 — g); and 3) for g < 1 has a cutoff at 
a; = — 1/(1 — g), below which it becomes identically zero. 

Using the function ( |T5i) . the basic expressions ( |T2l) and ( |T3l) for the OLM canonical 
formulation can be rewritten as 

p = Z-^e/(^-^^), (16) 

and 

= T,e-^(^-^'). (17) 
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It is then immediate to see that, for g — > 1, Eqs. (fT6|) and (fT7|l reproduce the conventional 
Boltzmann-Gibbs framework. 



III. TWO-TIME GREEN'S FUNCTIONS AND SPECTRAL DENSITY IN QUAN- 
TUM NONEXTENSIVE STATISTICAL MECHANICS 

A. Definitions and spectral properties 

In strict analogy with the extensive case , we define the two-time retarded (z/ = r) 

and advanced [v = a) q GFs in quantum nonextensive thermostatistics for two operators A 



and B as 



g^IbM = -reAt-t')mt).B{t%) 



1 



= mt)-B{t')))^^\ {u = r,a). (18) 

Here, Or {t - t') = 9 {t - t'), Oa {t - t') = -9 {f - t), and 9{x) is the step function. In Eq. 
f fTSj) . [..., ...]^ denotes a commutator (r/ = — 1) or an anticommutator {rj = +1) and X(t) = 
^iHtj^^-iHt -g ^YiQ Heisenberg representation of the operator X, satisfying the Heisenberg 



equation of motion (EM) 45 1 



^ = ^[H,X{t)]_. (19) 

The definition (ITSl) reproduces the conventional extensive formalism for q = 1 and allows 
us to develop the q GF framework equivalently with commutators or anticommutators. 
However, in practical calculations it will be convenient to use, in Eq. f llSp . rj = —1 or 
?7 = +1 for bosonic or fermionic operators, respectively. 

Physically relevant quantities, which enter the definition of q GFs, are the two time q 
correlation functions (GFs) FqAB{t,t') = {A{t)B{t'))q and FqBAit' ,t) = {B{t')A{t))q for the 
corresponding operators. Working within the equilibrium statistics, one can easily prove 
that the two-time q GFs and q GFs depend on times t and t' only through the difference 
T = t — t' . So one can write 

G%{t - t') = {{A{t - t')- = {{A- B{t' - t)))M. (20) 

This feature allows us to introduce the Fourier transforms 

^ ^G%{u;)e-^-\ (21) 

11 



2^ 



(22) 



where G 



qAB 



[uj) = {{A{t)] B))^q}u will be named the z/-type q GF of A and B in the u- 
representation and Fqxvi!^) = {^{.^)Y)q,iJj will be called the g-spectral intensity of the 
time-dependent q CF Fgxy(r), with ^{u) = /J^^ rfrS5(r)e*^" = G[%,Kar, . . 
We now define the time-dependent q SD for the operators A and B [2, y, [3] as 

A,ABir) = {[AiT),B]^),. 

For its Fourier transform (the q SD in the cj-representation) 

AgABiu;) = ([A(r),5]^),,,, 



(23) 



(24) 



one immediately finds the exact result 

. 2^ 



A,ab{u;) = {[A,B]) 



(25) 



which constitutes an important "sum rule" for the q SD AgAB(t^) to be used for physical 
consistency of practical calculations and approximations. 

Then, using the integral representation for the z/-step function 



+ 00 



dx 



0' 



(26) 



27TX + {-lyie' 

where the symbol (—1)*^ stands for +1 if z/ = r and — 1 if = a, we obtain the g-spectral 
representation 



for the Fourier transforms of the two-time q GFs (1201) . 

Combining Eq. flTTI) and the sum rule fl2S]) . one can get another general result which, as 
in the extensive case, may play a relevant role in practical calculations involving the q GFs. 
Indeed, as |a;| — > cx) we have 

du' AgAB{(^') 



-1 



27r 1 



([A 5] J 



rill 



and hence 



+ - / ^A.ab(^') K - (-l)'^^^] + O - , (28) 

UJ J-OD / 



uj-\ if ([AS],),^o 

(a>2), if ([A,5]^), = 0, 



(29) 
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which provides a relevant boundary condition for the q GFs. 

As in the quantum extensive case 0, [sl one can easily show that the q GFs G^^\^{uo)^ 
analytically continued in the complex uj plane, are analytical functions in the upper (for 
V = r) and lower (for v = a) half planes. Then, these functions can be combined to 
construct the q GF of complex u: 

which is analytical in the whole complex u plane with a cut along the real axis where 
singularities may occur. 

It is worth noting that, in terms of the SD, no formal differences exist for the spectral 
representations of the GFs in the extensive and nonextensive contexts. Hence, most of 
the developments in the extensive two-time GF framework remain formally valid in the 
nonextensive one. For instance, using the 5-function representation 

S(x) = lim —i— —\, (31) 

£^0+ 2m {x — 16 X + le j 

we have the relation 

AgAB{i^) = i[GqABiuj + ie)-GqAB{i^~'i£)]=i G^ab (^) " '^Sb (^) ^ (32) 

which expresses the q SD in terms of the related two-time q GFs in the u representation. 
This result, which is expected to play an important role in the applications of the q GF 
method (as happens in the extensive case), suggests also that the cut for GgAB^^^) along the 
real axis is determined by Eq. fl32|) and its singularities are given by real u values satisfying 
the condition AgAB^^^) 7^ 0. It is worth noting that, in general, GqAB^^^) has to be regarded 
as a many-valued function of the complex variable u. Hence, its singularities lie on the 
real axis on the first Riemann sheet. On the other sheets, the singularities may shift to the 
complex plane, leading to the appearance of complex poles. 

Another important result can be easily obtained assuming AqAB{^) real and using the 
relation 

lim du'- ^^""P . = V doo'^l^ - i-iytTrfioo), (33) 

where V denotes the main part of the integral. From Eq. fl271) . we obtain indeed 

duj' AqABji^' 

2n uj' — ui 



. ^ ^ l^"^ duj' AqABjuj') (-l)^ 

G\ab{^) = -V / iz-—, ■ —iAqAB{u}). (34) 
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Hence we get 

'qAB 



^^G),ab{^) = (-1) ^ / . (35) 



with 



A,^b(^) = -2(-l)'^ImG(';)^(^), (36) 



Kab{uj) = -2lmG'2B(^). (37) 



and, in particular, 

In analogy with the extensive counterparts , the relations fl35l) between the real and 
imaginary parts of G^^\q{uj) and ^^"^^(0;) will be called gdispersion relations or nonextensive 
Kramer-Kronig relations. 

In the next subsection we will derive spectral decompositions for KqABiy^), GgAB^^), 
FqAB{T) and FqAB{uj) which allow us to obtain information about the nature of the GF 
singularities and hence about the excitations in nonextensive quantum many-body systems. 



B. g-spectral decompositions 



Let {\n)} and {En} be the selected eigenvectors and eigenvalues of the Hamiltonian H of 
a many-body system and assume that {\n)} is a complete orthonormal set of states. In this 
representation, the q SD Ag^^l^), as given by the Fourier transform (|2^ . can be written as 



KqAB^uj) = 1^ J] [1 - (1 - g) /3 (E„ - [/,)] 

1 n,m 



q 

1-9 X 



X < 1 + ?7 



l-{l-q)(3{Em- Ug 



1 

1-q 



(38) 



[l-il-q)l3iEn-Ug) 
where use is made of the OLM canonical framework for calculation of the gaverages. In Eq. 

m, 

Z, = Trp'^ = J2[l-il-q)PiEn- f/,)]^ , (39) 

n 

^nm = {n\X\m), and Umn = Em — En- Besides, the spectral representation (!30|) for GqAsi^^) 
yields 

271 



n,m 



X <l + r] 



l-{l-q)P{Em-Uq 
l-il-q)l3iEn-Uq 



I -^nm.J-^n 



'■nm-^nm 



UO — UJrr 



(40) 
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Analogously, for the two-time g-CF FqAsij) and its Fourier transform FqABiyj)^ one easily 
finds 

FqAB{r) = S ^ [1 - (1 - g) /5 (E„ - f/,)]T^ A^^Bmne-''^-^ (41) 

n,m 

and 

FqABiuj) = S ^ [1 - (1 - g) /? (E„ - Ug)]^^ AnmBmnS {uJ - Umn) , (42) 
y n,m 

with similar expressions for FgBAiT) and FgBAi^^)- 

The comparison of the previous relations with the corresponding extensive ones 
indicates that the Tsallis statistics does not modify the meaning of the GF singularities, but 
changes substantially the structure of the spectral weights with the introduction of a mixing 
of the energy levels which is absent in the extensive framework. Equations (140|) and (14T|) 
suggest indeed that the real poles of GgAsi^^) (i-e., the frequencies ujmn) which are related to 
the eigenvalues of the Hamiltonian, represent the frequency (energy) spectrum of undamped 
excitations [oscillations in time of FgAsi'^)] in the system. It is worth noting that, for a 
macroscopic system, the 6 poles in the spectral representation (HOl) are expected to be lying 
infinitesimally close, therewith defining a continuous function Aab{^) for real uj. Then one 
can speculate that the excitation concept may work only under the basic assumption that the 
q SD exhibits some pronounced peaks whose widths have to be considered as a direct measure 
of the damping or the life-time of excitations (or gquasiparticles, elementary or collective 
depending on the physical nature of the operators A and B). As mentioned before, this 
picture should be associated with the appearance of further complicated singularities of the 
gGF GgAsi^) which may occur in the uj complex plane on the Riemann sheet below the 
real axis where G^^^q{uj) is not an analytical function. Hence, in practical calculations, one 
must search for the complex poles of G^^]^^{uj) very close to, but below, the real uj axis. For 
each of them, the real part will determine the frequency of the excitations (the q excitation 
dispersion relation) and the imaginary part will represent their damping or life time. In 
this scenario, AgAsi^^) will result a superposition of quasi- Lor entian peaks, at characteristic 
frequencies, whose widths will represent the damping of the related excitations. Of course, 
if these widths reduce to or are zero under appropriate physical conditions, the q SD will be 
given by a superposition of 6 functions signaling the occurrence of undamped excitations in 
the system. 
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C. Expressions of the two-time q correlation functions in terms of the g-spectral 
density 



In the previous section compact relations between G^^\^{uj) and KqAsi^) have been ob- 
tained in strict analogy with the quantum extensive counterpart 2|, [3| . We will show below 
that the peculiar nature of the Tsallis probability distribution prevents us from expressing 
the two-time q CFs {A{T)B)q and {BA{r))q directly in terms of the q SD KqABiyj) or the 
related q GFs G^q\^{uj) (z/ = r, a). 

First, it is worth recalling that in extensive quantum statistical mechanics it is a remark- 
able feature the existence of direct relations which allow us to express {A{t)B) = (A(r)i?)q=i 
and {BA{t)) = {BA{T))g=i in terms of Aab{^) or G^^(c<j) {u = r,a) for two arbitrary op- 
erators A and B. For future utility we remember that these relations read 

m 

and 

duo KAB{uj)e-''^^ 



iBAir)) = / — , (44) 



oo 



from which the corresponding static CFs can be immediately obtained setting r = 0. 

The situation becomes sensibly more complicated within the nonextensive context. To 
see this in a transparent way, it is convenient to start with the static q CFs for two arbitrary 
operators. Taking properly into account the presence of the 5 functions in the spectral 
decomposition (l38l) for the q SD KqABiyj)^ we can also write 

1 AgAB(t^) 



27rl I 



Veq 

n,m L + V^q J 

where we have conveniently introduced the modified q exponential function 

e:=[e^,V=[l + {l-q)x]^^. (46) 
Then, the integration over u easily gives the exact relation 

du AqABi^^) _ 



2vr 1 + r]eq 
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with Zq = TrCq 



By inspection of the right-hand side of this equation, one argues that it does not reduce 
to the correlation function {AB)q as happens in the extensive case q = 1- Nevertheless, if 
we introduce the "g operators" Ag and Bg (related to the original ones A and B) defined by 
the matrix elements 

Aqnm = {n\Aq\m) = Cq (n, iTi) Anm, (48a) 
Bqnm = {n\Bq\m) = Cq (m, n) Bnm, (48b) 



where 



with Cq=i (n, m) = 1, Eq. (gT]) yields 

m^r'£^. (50) 

J-oo 1 + r]eq 

Similarly, with AqBA^^) = V^qAsi—^), one finds 

iBM,^r^'^. (51) 

J-oo Cq + 1] 

The previous spectral relations, which express the q averages for products of the two q 
operators Aq and Bq in terms of the single q SD AqABi^^), are remarkably similar to the 
extensive static CPs for A and B and, as expected, they reduce consistently to them for 
g 1. Of course it is {XqYq)q ^ {XY)q (X, Y = A, B) ioi q ^ 1. 

Concerning the dynamical q CPs for q operators, using Eqs. fj42l) . fH5|) . and fH9|) it is now 
easy to show that the following relations are true 



duj AqAB{uj)e ~ 
2vr 1 + neq~ 



{Mr)Bq)q = I - T 1-,. (52) 



and 



{B,A,(r)),= I '^ ^'l^tr" ■ (53) 



oo 



which reduce to Eqs. (1431) and fl44p as g ^ 1. 

The substantial difference with respect to the case g = 1 lies in the unfortunate feature 
that, for g 7^ 1, the previous compact relations express the g averages for products of the 
complicated g operators Aq and Bq in terms of AqABiy^)- This may constitute a serious 
difficulty in exploring physical cases involving directly the CPs of the operators A and B 
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which enter the definitions of G^^J^ and AqAB{^)- Nevertheless, it is possible to obtain 
explicit relations, although cumbersome and in general not very handy, which relate the 
CFs for the physical operators A and B to those for the corresponding q operators. Indeed, 
from Eqs. (H71) . ( l50l) . and (ISB it is immediate to show that, for the static case, 

{AB),= rp-^^^^^-Z;^j:e-^^^^^^^^ (54) 
J-oo l + nen ~Z 



and 



{BA), = r - Z^'Y.~^,'^'"^-''''D,{n,m)B^^A^^^ (55) 

y ' / n.m 

with Dq{n,Tn) = Cg{n,m) — 1 — as g — 1. Similar expressions are true for dynamical 
CFs {A{T)B)q and {BA{T))q, which will involve exponential oscillations in time. 

A comparison with the corresponding extensive ones (H3|) and ( H^ allows us to under- 
stand, in a transparent way, the nature of the deviations from the extensive limit q = 1. It 
is worth noting that, under the condition 

|(l-g)/3(i/-f/,)| «1, (56) 

the previous relations simplify to 



-tWT 



and 



r^^f^ (57) 
(«^M),^r^^^4^, (58, 

J^oo Cq + 1] 

which become exact for q = 1. Of course, the simplified static q CFs for A and B can be 
obtained setting r = in Eqs. fl571) and fl58l) . These equations may be conveniently used in 
practical calculations. Of course, without the restrictive condition (j56l) . one must use the 
cumbersome Eqs. ( 15^ and ( l55i) and resort to numerical calculations taking properly into 
account the cutoff condition ^\M . It is worth noting that, although the condition ( |56l) is 
certainly verified in the limit g ^ 1, it can be also realized for q far from unity with suitable 
choices of the parameters (3 and En — Uq. However, in view of a still reduced number 



of applications |27l . l46l |. in order to gain experience in using the q many body formalism 
developed before for more complex situations, it may be in any case useful to consider weak 
nonextensivity conditions with q close to unity. 
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D. Parametric representation for the two-time q Green's functions 



As a conclusion of this section, we will present an alternative way to introduce the two- 
time GFs, the related SD , and the two-time CFs by using a parametric representation, 
suggested in Refs. 



25 



26[ | in the context of the Kadanoff-Baym formalism This rep- 



resentation allows to express the relevant q quantities in terms of appropriate parametric 
integrals involving the corresponding q = 1-ones. 

The aim is to clarify the statement given in the Introduction about the effectiveness of 
the method in practical calculations and to stress again the potentiality of our framework 
for a direct calculation of the gquantities of interest. The basic idea is to take 6=1 — 
{1 — q) j3 {H — Uq) and alternatively z = 1 + 1/(1 — g), z = 1/(1 — g), and z = q/ {1 — q) in 
the contour integral representation 

du 



ih 



l-z 



2n 



exp(— m6)(— n) 



(59) 



with 6 > and Re^; > 0. Here C denotes the contour in the z complex plane which starts 
form +00 on the real axis, encircles the origin once counterclockwise, and returns to +oo. 
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26| the following representation for ^q(/3), 



With these ingredients, one can easily obtain 
^'qAB i'^'i f^)i ^^"i {^q{'^)Bq)q (hcrc wc uccd to explicit the /^dependence) 



UP) 



duK^^\u)Z^ [-(3u{l-q)], 



GtLij; P)= [ duKf\u)Z, [-Pu (1 - q)] G[% (r; -(3 
Jc 



and 



where 



u{l-q)) 



{Ag{r)Bg)g= / rfwi^f (^i)Zi[-/3n(l-g)](A(r)i?)i,_^„(i_,), 



c 



u 



-T 



271 

Kf\u) = 
Kf\u) 



l-q 



1 \ g-«[l + (l-q)/3C/,] 



-u) ^-1, 



l-q 



u 



(60) 
(61) 

(62) 

(63) 
(64) 
(65) 



and Zi, G^^}^, and {A{T)B)i^j3 denote the corresponding extensive quantities. 

An analogous integral representation for the q SD A^ab (^^; P) can be simply obtained by 
replacing, in the Fourier transform of Eq. (IMl) . the relation (^7^ which connects 0^!^^^ {u; (3) 
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to KqAB (i^; /3) and, for g = 1, G\'Iq {oj; [3) to tlie extensive spectral density Ai^^ (cu; [3) = 



In view of tlie previous cumbersome (altfiougli elegant) parametric representation, in our 
opinion the method to reduce q many-body problems to the corresponding extensive ones 
does not appear, in general, convenient in practical calculations. The two main reasons 
are (a) it involves the a priori explicit calculation of {q = l)-quantities which, except for a 
limited number of simple cases, requires a first step of more or less reliable approximations; 
(b) after that, one must calculate nontrivial contour integrals of complicated functions. In 
general, this may be a formidable problem which could require additional approximations. 

Thus, we consider it worthy of interest to develop many-body methods for a direct study 
of the q properties overcoming the a priori knowledge of the related extensive ones. We 
introduce below the appropriate extensions of two well-known and powerful quantum many- 
body methods (for the classical framework, see Appendix A). Next, we will support our 
statement by means of a direct detailed study of the g-induced nonextensivity effects on the 
low-temperature properties of a nontrivial many-boson model. A preliminary study, along 
this direction, of a c?-dimensional Heisenberg spin model with long-range interactions was 
performed in Ref. j^]- These are only two nontrivial many-body problems used by us to 
test the effectiveness of our suggestion to perform direct nonextensivity calculations within 
a genuine q many-body theory. 

IV. METHODS FOR DIRECT CALCULATION OF THE TWO-TIME q GREEN'S 
FUNCTIONS AND THE q SPECTRAL DENSITY 

In this section we present an extension to the nonextensive quantum statistical mechanics 
of two intrinsically nonperturbative methods in strict analogy to the extensive counterpart. 
These methods will be called the q EMM, for a direct calculation of the two-time q GFs, 
and the q SDM, for a direct calculation of the g-SD. In principle, in view of the exact 
relations established in the previous section, both the methods are completely equivalent 
in the sense that they should give exactly the q GFs and the related q SD. Nevertheless, 
previous experiences in quantum [7|] and classical |9| extensive statistical mechanics suggest 




. One finds 



A,AB (uj; (3) = I dwirf Zi[-/3M(l-g)]AiAB(a;;-/3M(l-g)) 



Jc 



(66) 
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that, in practical calculations, the g-SDM may have several advantages for making more 
systematic and controllable approximations. 

A. q equations-of-motion method 

By differentiating Eq. (!20|) with respect t = t — t' we obtain the q EM for the q GFs in 
the r representation, 

{{A{r)- B))^^^ = 6 (r) {[A, B]^) ^ + {{[A{t), H]_ ; B))['^ , (67) 

or in the a; representation (more convenient in practical calculations), 

^ ((^(r); = {[A, B]^)^ + {{[A{t), H]_ ■ B))^^ . (68) 

Equations ( l67l) and ( l68ll are not closed because higher-order q GFs occur in the problem. 
Therefore, one needs to consider an additional EM for these new functions which is again 
not closed. By iteration of this procedure, we obtain an infinite hierarchy of coupled EMs 
of increasing order which can be written in a compact form as 

{{L-Air); = 6 (r) {[L^A, B]^) ^ + ((^^^^(r); B))[''' , (m = 0, 1, 2, ...), (69) 

or 

iv {{L-Air); B))'^l = {[L^A, B]^) ^ + {{LT'A{t)- B))^^1 , (m = 0, 1, 2, ...), (70) 

in time and Fourier representation, respectively. Here, we have used the EM f fT9|) in the 
form 

^^^A{r) = L'^A{T) = [A{r\H]_, (71) 

and the operator means L^^A = A, L\jA = [A, H]_, L^A = [[A, H]_ , H~\_, and so on. 
Note that the chain of q EMs in the representation ( |69l) or (!70l) is formally the same for 
different types of q GFs and hence one can eliminate the index u when the physical context 
does not offer ambiguity. 

To solve the chain of EMs in the form (169|1 and (ITOjl . we must add appropriate boundary 
conditions which, in the cjrepresentation, can be identified with the asymptotic behaviors 
(j29l) . Of course, although Eqs. (!69l) and (jTOll are exact, it is impossible to find a complete 
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solution for interacting systems. In practical calculations one is forced to use decoupling pro- 
cedures, and hence approximate methods, to reduce the infinite chain of coupled equations 
to a finite closed one which may be solved. However, in general, systematic and controllable 
decouplings are not easy to find and one must check the reliability of a given approximation 
by comparing the results with experiments, simulations or other types of approaches. The 
gSDM, which will be the subject of the next subsection, should be more flexible in such 
direction as it happens in the extensive case 



B. (/-spectral density method 

By successive derivatives of Is^qAsij) [Eq. (!23|) ] with respect to r and using the EM (ITTIl . 
we have 

^A,AB(r) = Hr([LSA(r),i?],)^ (m = 0, 1, 2, ...). (72) 

Then, taking the Fourier transform of Eq. ( |72i) and setting r = 0, integration over uj yields 
finally the infinite set of exact equations for KqAsi'^)'- 

/ + 00 J 

= {[A,C^B],)^, (m = 0,1,2,...), (73) 

where the operator means C%B = B, C\jB = [H,B]_, CjjB = [H, [H,B]_'\_, and so 
on. The quantity on the left-hand side of Eq. (!73l) will be called the m moment of AgAsi^), 
and the relations ( ff3l) constitute an infinite set of exact moment equations (MEs) or sum 
rules for the q SD. 

The infinite set f l73l) can be seen in a different way. Since the rj commutators and hence the 
q expectation values involved on the right-hand side can be calculated, at least in principle, 
it is quite remarkable that the m moments of the q SD can be explicitly obtained without a 
priori knowledge of the function AqAB{^)- This important result implies that the sequence 
of Eq. ( |73|l represents a typical moment problem. Its solution would yield the unknown q SD 
and hence all the related quantities (g GFs, q CPs and other observables) . Unfortunately, 
also this problem cannot be solved exactly and one must look for approximate solutions 
along the lines specified below which constitute the key idea of the original SDM Q]. 
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1. Polar ansatz 



As suggested by the exact spectral decomposition (1381) . one seeks for an approximation 
of AgAsi^) as a finite sum of properly weighted 5-functions of the form {polar ansatz) 

n 

A,M^) =2nY, - J^), (74) 

where n is an integer number. The unknown parameters A^*^^ and t^g^^^, depending on the 

physical nature of the operators A and B, have to be determined as a solution of the finite 

set of 2n equations obtained by inserting expression (1741) into the first 2n MEs, Eq. (1731) . 

(fe) 

Physically, the parameters oj^aq play the role of effective eigenvalues of the Hamiltonian 
and each of them represents a real pole of GqAsi,^) corresponding to a mode of undamped 
oscillations for the gCF {A{t)B)^ [see Eq. (j4T|)]. 



2. Modified Gaussian ansatz 

As outlined at the end of Sec. IIIIBt there are physical situations where the damping of 
oscillations in the system under study may be relevant and hence the polar approximation 
f l7^ is inadequate. In these cases, the basic idea of the SDM, related to the moment problem 
( I73|) . remains still valid, but it is necessary to choose a more appropriate functional structure 
for the q SD which allows one to determine the modes of excitations in the system and their 
damping or life time. In the extensive context, a generalization of the SDM in this sense 



was first proposed by Nolting and Oles [47| for Fermi systems and by Campana et al. [48] 
for Bose and classical systems whose SDs are not positive definite in the whole range of ui. 
Following their key idea, to assure the convergency of the q SD moments at any order and 
to preserve the intrinsic physical character of KqAsi'^)-, one can assume for the q SD the 



modified Gaussian ansatz 



48| 



A,ABiu;) = 2n{ef- + r,)J2^^e ^'^s . (75) 

Clearly, with the functional representation fl75l) for AgAsi^^), the width of the peak in 0; = 
'^qAB related to the parameter F^^^ and the life-time of the excitations with frequency 



ujgAB to be identified with t^'^^ = y^l^AB under the condition T^^ab/ 
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(fc) 

'^qAB 



2 



< 1. The 



choice ( !75|) is only motivated by the fact that it makes direct contact with the notation used 



in the hterature for extensive problems 



47 



48l | and, in view of previous experiences, it is 



expected to simplify the algebra also in explicit calculations about g-induced nonextensivity 
effects. 

As in the q EMM, also in the q SDM the problem remains to close the truncated finite 
set of q MEs arising from the polar ansatz (17^ or the modified Gaussian ansatz (!75|) . In 
any case, evaluation of the right-hand side of Eq. fl73l) should generally involve higher-order 
q SDs. Hence, higher-order moment problems should be considered and the difficulty of 
calculations will increase considerably. So, in order to solve self-consistently the finite set 
of q MEs, which arises from Eq. (175]) using the ansatz (1711) or (175]) . it is usually necessary 
to use some decoupling procedures and thus to introduce, in a systematic way, additional 



consistent approximations in the SDM as in the extensive case 



fly, 



48|. 



V. NONEXTENSIVITY EFFECTS FOR A HIGH-DENSITY BOSE GAS WITH 
STRONG ATTRACTION BETWEEN PARTICLES 

A. Model 

For practical and explicit calculations we consider here a nontrivial Bose model introduced 
several years ago by Babichenko 3l| with the aim to explore the properties of a Bose system 
with strong attraction between the particles. In contrast with real situations, where the 
interparticle interaction potential is characterized by a hard repulsive core at short distances 
and a strong attractive well at large distances, in Ref. jsil the repulsive core was assumed 
to be soft. This simplification allowed us to study exactly the properties of the model at 
T = 0, making it possible to use diagrammatic techniques quite similar to the ones used 



for high-density Bose systems with a Coulomb pair interaction [49|]. It is worth noting that, 
although the choice of a soft repulsive core does not correspond to the real situation, the 
possibility to obtain exact results for such a model is undoubtedly of interest since it leads 
to a strongly compressed ground state strictly related to the peculiar relation between the 
parameters of the attractive and repulsive parts of the pair interaction potential. In this 
sense, the model may be considered as a complement to the well-studied weakly nonideal 
Bose gas. 
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The same high-density Bose model was further studied in Ref. [32] to explore the finite- 
temperature effects in the context of the extensive many-body theory. Specifically, the usual 
Bogoliubov approximation was proved to be valid for this model and the exact results by 



Babichenko |3]| were simply reproduced. Besides, the low-temperature properties of the 
model were again achieved without employing cumbersome diagrammatic techniques. 

On this ground, also as a further contribution to the general Bose-Einstein condensation 
(BEG) scenario, we apply the g- many-body methods of Sec. [IV]to investigate the g-induced 
nonextensivity effects on the low-temperature behavior of the relevant thermodynamic quan- 



tities for the Babichenko model 3l[ |. 



As a first step, we present below the definition of the model. Working in the grand- 
canonical ensemble and with periodic boundary conditions, the Bose model of interest is 
described by the second- quantized Hamiltonian 



n = H-fiN= / (frij\r) I — ^ - /i \ ^(r) + 



2 

where N is the total number operator of spinless bosons with mass m, /i is the chemical 
potential, and '?/'(r)and ip^r) are the usual Bose field operators. The pair interaction poten- 



tial ip (r) in fl76|) is assumed [3l|] as the superposition of a repulsive Yukawa-like potential of 
radius Rq and of an attractive Gaussian well of radius R > Rq and depth Uo > 0, with the 
representation 

^(r) = ^exp(-^)-f/oexp(-^). (77) 

Here 7 is a certain definite positive coupling parameter and the strong attraction or deep- well 
condition UqR^ > ff/m is assumed to be satisfied. 

For our purposes it is convenient to choose the system of units with H = m = Rq = 1 and 
to work in the wave-vector {k} representation. So the grand-canonical Hamiltonian (!76|) 
assumes the form 

^ = Xl^k^ak + ^ (|ki - kal) (5ki+k2,k3+k4 4i«k2«k3ak4, (78) 

k {k,} 

where = — fi, V is the volume of the system and the Fourier transform ip (k) of the 
pair interaction potential (!77j) is given by 

^if^) = TT¥- -'^'ua^ (-^) ' (79) 
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with 

(0) = 7 - n^/'^UoR^ = 7^=^/2r/f/oi?^ < r/ < 1, (80) 

assumed to be small. 

With the previous definitions, the coupling parameters involved in the pair interaction 
potential are connected by the relation 

^ = {l + r])iT''/^UoR^ (81) 

and from the strong attraction condition UqR^ ^X±,one easily finds that 7 > 1. 



It is worth mentioning that, as shown in Ref. [SI] and confirmed by the following calcu- 
lations, the Gaussian form of the attractive part of (f (r) is not essential. Indeed, the key 
condition to be used through the calculations for the second term in Eq. flTTl) is only a 
sufficiently smooth change with distance of (r) such that its Fourier transform is localized 
in a small region of the wave- vector space {k < 1/R)- 



B. q equations- of- mot ion method within the Bogoliubov approximation 

Adopting the conventional Bogoliubov approximation, usually restricted to weak inter- 
actions and low densities [s^, but proved to be valid also for high-density charged Bose gas 
5l| and for the Babichenko high-density Bose model with strong attraction between the 



particles 



32| . the grand-canonical Hamiltonian (!78|) reduces to 



1 N'^ 



(82) 



where NgQ denotes an unknown q mean number of bosons in the condensate and 



fgk = + ngo[Lp{0) + (p (k)], 



k 



rigo^ (k) , n, 



N, 



qO 



qO 



V 



(83) 
(84) 



According to the standard Bogoliubov picture, one now should diagonalize the truncated 
Hamiltonian ( l82l) by a linear canonical transformation and then proceed, at least in principle, 
to calculate the OLM low-temperature g-thermodynamic quantities with the prescription 
outlined in Sec. [Ill (for the extensive case see Ref. (32|). However, we find it convenient 
to follow a different approach which makes direct contact with the q many-body formalism 



26 



developed in Sees. IllltHVl and allows us to explore, via direct explicit calculations, the q- 
induced nonextensivity in the low-temperature regime in a simpler and more transparent 
way. 

Let us introduce the single-particle retarded and advanced two-time q GFs 



(85) 
(86) 



where the Heisenberg representation of the operators and the q averages have to be consid- 
ered with respect to the truncated Hamiltonian fl82|) . 



Working in oj space, for the retarded or advanced q GFs Gqy_ {oj) 



ak(r);ai 



and 



Ggk (t^) 



one easily finds the two coupled EMs 



t^Ggk {Uj) = 1 + fqkGqk (Uj) + hgkGqk (uj) 
LoGqk (t^) = -hqkGqk {uj) - fq^Ggk {u) . 

This algebraic system can be simply solved to find 



Gqk (t^) 



1 



/gk 



1 



^qk J ^ — ^q\i 



1 



^qk / + t^gk 



Gqk {uj) = 

Here, the quantity 

t^gk = (/gk ~ ^gk) 



igk 



gk 



- ^gk 



2ujq]i \UJ — tUqk to + tUqk 



(87a) 
(87b) 

(88) 
(89) 



H 2 



- fJ' + nqo {ip{0) + '^{k)) 



-^gO<^^(^) 



(90) 



represents the energy (frequency) spectrum of the undamped elementary excitations in the 
system. To obtain the relevant q physical quantities it is now convenient to introduce the 
single-particle q SDs for the two gGFs Gq^ and Gq^ defined as [see Eqs. ( I23l) and (E 



Agk(t^) 



ak(r),ai 



and 



respectively. 



Agk(c^) = ( a\{T),a 



q,uj 



(91) 



(92) 
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Then, from Eqs. (l32l) . (l88ll . and (l89ll . we immediately have the two (5-function represen- 
tations for the spectral densities: 



Agk(c^) = vr [(1 + 7gk)5 (t^ - u;gk) + (1 - lqk)S + cUgk)] 



(93) 



and 



where 



fqk , /i, 



7gk 



A 



gk 



gk 



(94) 



(95) 



We are now in the position to formally determine the relevant g-thermodynamic quantities 
simply taking into account the general relations stated in Sec. Illll 

First, we focus on the chemical potential that, within the Bogoliubov scenario, is defined 
by I52I /i = ZdH/dNgo) . Since H = n + fiN, with iV ~ AT^q + Ek^o "k^k, Eq. ([82]) yields 



kT^O k^O L 



(96) 



Then, with simple algebra, we get 

kT^O 



Besides, for the q internal energy Ug = (H) , we find 



(/gk - ^k) (aiak) + (ctk'^-k 



(97) 



k^O 



/gk ( aJc'^k ) + /igk ( aj^a 



(98) 



In Eq. ( 198!) . = (^^^^ and — Ai'go = J2kj^o (^'^k'^k^ , where is the total number of 
bosons in the system. In particular, the q depletion of condensate will be given by 

^E(4«k) , (99) 



n — n, 



gO 



k^O 



where n = N/V is the total number density of particles. Notice that, using Eq. fl99l) . 
expressions ( 1971) and (1981) for /x and (h\ can be also written in the most compact form 



k^O 



aj^ttk 



(100) 
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and 

(n) = iV,o(/. - n,o^(0)) + i^S^ + 1 ^ A:^ (4,,) . (101) 

Due to the complicated relations and (155|) . it is not easy, in general, to calculate 
explicitly the q averages which enter Eqs. fl971) - fl99l) in terms of the q SDs Agk(ci;) and 
Agk('^). Intricate numerical calculations become necessary, or one must resort to reliable 
approximations to obtain analytical results. The problem becomes sensibly handier under 
the condition [see Eq. 

'[l-q)l3(H-U,) <1, (102) 



with Uq = Ug — fi (^^^ • Under this condition, from Eq. S57}i it immediately follows that 

\ / q J -00 ^TT Ca - 1 

f + OO 



9 

and _ 



On the other hand, in the present Bogoliubov scenario, kq^jyoj) and l^qvioS) are exactly 
expressed by Eqs. ( 193|) and ( 19^ . Then, from Eqs. (11031) and (I104p we easily find 

-W) = - 1 + I (105) 



and 



where and Agk are given by Eqs. (p5|) . 

To calculate explicitly and consistently the g-thermodynamic quantities as functions of the 
temperature T and the total number density n, one must now express the chemical potential 
as a function of T and n, performing the sums over k in the previous basic expressions. This 
is not an easy task since, as q, fi enters the problem in a cumbersome way through the energy 
spectrum cj^k- However, it will be shown below that, in the physical regimes of interest (high- 
density and low-temperature regimes), one has |/i — /Xgol //^go ^ with figQ = UgQip^O), and, 
consistently, (n — Ugo) /n = [N — Ngo) /N <^ 1, as expected in a Bogoliubov approximation 
scenario. This key feature allows us to simplify sensibly the problem setting, as a good 
(first) approximation, /i ~ /igo = UgQip^O) = {ngQ'~^)rj / {1 + rj) and also UgQ ^ n in the 
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expressions of ^k, ^qk, and /^k, which enter the previous general relations for the relevant 
g-thermo dynamic quantities. So, in particular, with fi ~ fiqo + A/i, the correction Afi to the 
chemical potential will be given by the last two terms of Eq. (1971) with fi replaced by n!f{0). 
In all cases, for explicit calculations, the expressions of e^, cj^k and /^k in the summands 
can be approximated, to leading order, as 



/<?k 



-n(^(0), 
— + nfcV(k) 



+ rnp{k) = /k. 



(107) 
(108) 
(109) 



Notice that c<Jk, on the right-hand side of expression (llOSp . is just the Bogoliubov energy 
spectrum of elementary excitations. 

After that, with Ek^o(-) {271^)-^ J^^ dkk\...), very complicated integrals re- 
main again to be performed due to the presence of the power-law function involved in 
the integrands. Since resorting to numerical calculations by variation of (3 and q is beyond 
the purposes of the present work, to obtain demonstrative explicit results one is then com- 
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pelled to make further suitable approximations as the known factorization procedure 
or expansions close to the extensive value g = 1 of the Tsallis parameter. We follow here 
the second direction because it is simple and instructive to estimate, to the leading order 



in g — 1, the q nonextensive effects on the already known exact extensive results [3^ in the 
high-density and low-temperature regimes where the integrals can be explicitly calculated. 
To the first order in g — 1 one obtains the expansion 



1 



1 



eg-'- - - 1 e^'-"" - ± z (^e^ 

So, for the basic g-thermo dynamic quantities, we can formally write (with a cumbersome 
but simple algebra) 



;iio) 



n - UqQ 



2V ^ 



-E 



/k 



V ^ — ' uu e 



+ 



E 



y ^ (e'^-" - 1)' 



'nil 
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/i ~ rnp{0) 



1 

n 



1 

2n 



V ^ e^^^ - 1 



2 /k — t^k 



kT^O 



C^k 



+ 



+ 



A^Vk 



kT^O 



1 



+ 



;ii2) 



^\ 1 1 



t^k 



k^O 



kT^O 



(113) 

I i^l-'^k — II I 

k^O k^O 

Now all the integrals in Eqs. (11 lip - ( 11X31) can be consistently calculated under the con- 

1 /2 

ditions (wy) ^ 1 (high-density regime) and /5 ^ 1 (low-temperature regime) where it 
is reasonable to speculate that the Bogoliubov approximation preserves its validity 32|. 
Indeed, making the transformation x = k^wjY^'^, one sees that the principal contribution 
to the integrals is made by a; ~ 1 so that neglecting 0((n7) ' ) or exponentially small 
terms is quite legitimate 3l(|. A similar situation arises also in a high-density Bose gas with 



Coulomb interaction between bosons 



49, 



5l| . Bearing this in mind, in the high-density 



and low-temperature limits, for the depletion of condensate, the chemical potential, and the 
q internal energy density Ug = Ug/V = (^H^ /V as functions of T and n, one finds (see 
Appendix B for calculations of the integrals) 

'ri + V 



n - ngo ^ a (717) * + — 



(n7)"5 



+ (1-9) 



1 /r/ + l 



6 



r7 + 1 



TT 



(^7) -67(^7)4 + 

r] 60 



+ (1 



^' 15 



ri 

Tj + l 

V 

V+ 1 
V 



27r2 fr] + l 



{n^r'T}, (114) 



lr/ + 1 1 , ,2 4, , ,5 vr^ 
2 7 5 30 



7(^7) 
7(n7)"5r3 

. 3 

ri+1 



+ (1 



^^15 



v + i 



(n7) 



1 + 2 
1 + 2 



{n^r'T}, (115) 



r7 + 1 



(n7)-'T , (116) 
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where a = V^T'^{l/4)/{487i^/^) and b = V^T'^iS/A)/ {An^/^). 

Notice that, from Eqs. flll4p and (11161) . for (ngo7)^^^ — (^7)^^^ 1 and T ^ 1 one 
has (n — Ugo) /n and |/i — yU^ol //UgO ^ 1; as expected, signahng the internal consist ency 



of calculations. Of course, for q = 1 the extensive results are exactly reproduced 



31 



32|. 



From the previous basic relations one can now calculate all the g-thermo dynamic quantities 



of the system using the OLM formalism 
volume is given by 
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30l | . For instance, the q specific heat at constant 



3 



+ (1 - ^ V (-7)- T^ll + U'-^) in,r T } . (117) 



5\?7/ 3 \ rj 

It is worth noting that, as suggested by previous results, the g- induced thermal effects may 
compete sensibly with the extensive contributions. 



C. g-spectral density method at work: The two-pole approximation 

In this section we apply the q SDM directly to the grand-canonical Hamiltonian flTSl) by 
avoiding the Bogoliubov approximation and hence the known troubles related to a truncated 
Hamiltonian which no longer conserves the particle number and to the unwanted dependence 
of N,o. 

Taking advantage from the previous results, we will work rather via a two-(5-function 
ansatz for the appropriate q SD as suggested by Eq. fl93l) . Then, we will establish the 
conditions for which the results of Sec. IV Bl can be reproduced in the high-density and 
low-temperature regimes. As we shall see, the procedure will also open a window towards a 
possible systematic study of the g-thermodynamics of other more realistic second-quantized 
many-boson models. 

For the Hamiltonian (1781) in the wave-vector representation, we introduce the single- 
particle q SD 

A,kM = ([ak(r),4]_\ . (118) 
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Then, the general q MEs, Eq. (173|) . for the present problem can be written as 



[ak,£^ai]-) , (m = 0,l,2,..). (119) 



Using the usual bosonic canonical commutation relations, tedious but simple algebra yields 
for m = 0, 1, 2 (the case of interest for us) 



hoo 



oo 



duj . , ^ , ^ 

— A,k(cc^) = l, (120) 



/ duo 1 
— c.A,k {uj) = e^ + + - J] V'dk - k'|)iV,k', (121) 

— o^^Agk (^) = -el + 2£k y (^) + + ^ '^d'^ - k'DAgk', (122) 



where 



n 



k 

iV,k = (aiak)^, (124) 
Lq = E5ki+k2;k3+k4^^ (|ki - kal) (^aj^^flka^a^ki)^, (125) 



and 
Aqk 



^^¥^i\^-'^'\)Nqk' + V'dki -k|)(5ki+k2;k+k4 (4i42akak4) • (126) 

k' {ki,k2,k4} 



As we see, Lg and A^k introduce in the problem two-particle q CPs which should be expressed 
in terms of A^k (^) to close the truncated system of MEs fll20p - (11221) . Unfortunately, this 
constitutes a serious difficulty since, as shown in Sec. IIIIl there is no simple relation between 
a CP of the type {BA)^ and the corresponding q SD KqAB (^^)- However, under the condition 



'i-q)(3(n-u. 



<^ 1, the calculations simplify sensibly. Indeed, in this case, one can 
assume {BA)^ ^ (27r)"^ J^^ duKgAB (^) /(e/"" - 1) and in Eqs. ffT2T]) -ffT2D and ffM one 



Besides, with standard straightforward algebra 



7 



32| . we have also 



^ XI '^dki - k|)(5ki+k2;k+k4 (4^42«kak4)^ = j -^ -pcu ^''^ V (^) • (128) 



{ki,k2,k4} ^ ^-oo 
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This allows us to express Agk in Eq. (11260 in terms of Agk (uj) as 

"^'^ diu A„k, (iu) /•+°° diu cu-e^ 



and to obtain for the q internal energy the following expression in terms of the single-particle 
q SD: 



^ 1^ — c 



duj UJ + Ek 

^ ~fs. / .k(^), (130) 



where = X^k^gk- 

Nevertheless, the quantity cannot be expressed exactly in terms of Agk (uj), so that, the 
truncated system (I120p - (ll22p is not yet closed and one must resort to additional decoupling 



procedures 



a 



which allow one to express also Lq in terms of Agk (uj). We shall see below 



that this problem can be simply solved on 



hvsical grounds in the regime of interest. 



Working within the spirit of the SDM [7|, (see Sec. Ill), we assume for A^k (uj) the 



two-5-function ansatz [also suggested by Eq. fl93l) ] 

Agk(t^) = TT [(1 + 7gk)(5 {UJ - t^gk) + (1 - 7gk)5 (t^ + t^gk)] • (131) 

Then, with some algebra, Eqs. fll20p - (11221) yield, for the unknown functional parameters 
7gk and cUgk, the self-consistent equations 

7gkt^gk = £k + V{0)n + i ^ ,^(|k - k'DAT.k', (132) 

k' 

2 , , „ 2 



'^Jk = L, + T7 5Zv'(|k-k'|)A,k'+£k(£k + 2n^(0)) + -5^<^(|k-k'|)Ar,k. 
k' k' 



In these equations A^^gk is formally given by Eq. (11050 and 



(133) 



k' ^ 



1 J i^^qk - £k) (1 + 7gk) (t^gk + £k) (1 - 7gk) 



Cg J- Og ± 



(134) 



Of course, when 7gk and cjgk are known, for the q GF related to Agk (tu) 

+°°cia;'Agk (cu') 



Ggk(.)-/ (135) 



oo 
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iri + 7,k , l-7,k\ _ ^ + ^k + ri<^(0) + ^Ek'^(|k-k'|)Ar,k' 



we will have 

As expected, Eq. f ll36p implies two poles uj = icj^k (with cj^k > ) on the u real axis 
and hence the parameter ujgk determines the genergy spectrum of the undamped elementary 
excitations in the system, formally given by 



UJ 




-Jk=U, + |j]v'(|k-k'|)A,k' + ek(^k + 2n^(0) + |,J]^(|k-k'|)iV,kO!> . (137) 



^ k' 



In the polar ansatz fllSip . to determine Aq]^{uj) and hence all the relevant g-thermodynamic 
quantities, the problem remains to express the unknown higher order quantity Lq in terms 
of the parameters 7gk and cj^k- 

This difficulty can be easily overcome if we limit ourselves to explore the condensate 
phase (which is of main interest for our Bose model in the high-density and low-T limits) 
when a macroscopic population of the zero-moment single-particle state takes place. Indeed, 
in this situation one can obtain a formally exact expression for Lq setting co'gk=o = below 
a certain critical temperature , to find, from Eq. ( I137p . 



+ ( 2y;(0)n+ I 5^y.(A;')iV,k') - | V^(A:')A,k', (138) 
\ k' / k' 



where A^^^k and A^k are given by Eqs. fll27p and (11290 . 

Then, in the condensate state the q energy spectrum of the elementary excitations can 
be written in the form 

1 



where 




^5^y.(|k-k'|)iV,k' + -(/i-n^(0)) 



V 

k' 



+ fi,k^ , (139) 



= I E (|k - k'l) - (|k'|)] (A,k' - /iAT^k') . (140) 
k' 

Of course, Eq. (11391) does not provide the explicit g-dispersion relation as a function of T 
and n since A^^^k and A^k contain cjgk itself and 7gk which have to be still determined as 
solutions of the closed system of equations (132) and (I139p . 

In general, this problem is difficult to solve and, as for the extensive counterpart [tI. IssI. 
one is forced to consider asymptotic regimes for obtaining explicit results or to use numerical 
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calculations. Obviously, solving numerically the previous closed self-consistent set of coupled 
q MEs, by variation of T and q, constitutes a formidable tour de force (beyond the purposes 
of the present paper) which requires a separate study. Here, to avoid obscuring complicated 
calculations, we follow the procedure already used for the extensive case to obtain, in a 
natural way, predictions for the condensate state assuming 



N -N, 



qO 



n - riqo 



N 



<1, 



n 



and 



|A/i| 



< 1, 



(141) 
(142) 



where A/x = — rnp{0). 

These inequalities, already used in the previous subsection just as the basic conditions 
for the validity of the Bogoliubov approximation, must be consistently satisfied at the end 
of calculations. 

It is easy to check that under conditions (11411) and fll42p . with ~ nip{0) and n^o — n 



to leading order, the q MEs simplify to 



56] 



Iqk^^qk ~ k'^/2 + n(f{k) 

2 



LU 



qk 



Here 



(B) 



+ nA;V(k) 



(143) 



(144) 



is the Bogoliubov spectrum, and 



^<1^\ ^l~nip{Q) 



n 
V 



J2v{k') [^{\k-k'\)-^{k')]{i + 



(S) 
k' 



- 1 



(S) 

k' 



- 1 



+ 



1 C(jr 

provides an estimate of the correction to the Bogoliubov spectrum where the inequality 



(145) 



{B) 



1 2 



has to be checked at the end of calculations. Of course, from Eq. (11431) . 



we have for the parameter 



7<zk 



{e/2 + nip{k)) 



(B) 



(146) 
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Thus, the original g-moment problem is solved and all the g-thermodynamic properties for 
the condensate region under the conditions flMip and (11421) . and in particular the q depletion 
of the condensate n — n^o = XlkT^o ^9k, follow immediately from the spectral relations (11271) . 
(11301) . and (11351) . with Agk(ti;) given by the two-5-function representation (11311) . All the 
relevant expressions in terms of cjgk and 7gk at the Bogoliubov approximation level can be 
easily derived or also obtained from Ref. 55| with replaced by . These are quite 
cumbersome and not particularly instructive and hence they will not be reported here for 
brevity reasons. 

It is remarkable that, by a straightforward extension to the case g 7^ 1 of the procedure 
used in Ref. {ssl for the extensive case, one can systematically estimate, at least in principle, 
the corrections to the predictions of the Bogoliubov approximation for any plausible potential 

For the high-density Bose model under study in the regime of interest, the parameters 
cjgk and 7gk in all the summands can be replaced by their leading expressions uj\^ and 
7k^^ = [A;^/2 + nip{k)]/ uj^'^ and the basic equations become essentially identical to that 
obtained in Sec. IV Bi In particular, with (727) > 1 and /? > 1, all the results (|TTill -( |TT7j) 
near to the extensive regime are easily reproduced and the conditions (I14ip and (I142p are 
found to be consistently satisfied. 

As a conclusion, we make some comments which may be of practical interest for future 
calculations. For the high-density Bose model defined in Sec. IV Al and possibly for other 
models, the q SDM offers some advantages with respect to the method used in Sec. IVBI based 

on the Bogoliubov scenario. Adopting the q SDM one avoids the a priori replacement aj, 

1/2 

an — » A'^Q in the original Hamiltonian and hence some consequent conceptual difficulties 



57l | in both the extensive and nonextensive cases. Besides, it allows one to obtain, in a 



systematic way, the corrections to the Bogoliubov predictions. Finally, the procedure may 
be adapted to other more realistic situations provided that the conditions ( ]14ip and ( 1142^ 
are consistently satisfied at the end of calculations. 



VI. CONCLUDING REMARKS 

In this paper we have extended the general formalism of two-time GFs {2, 3, 7| to nonex- 
tensive quantum statistical mechanics in the OLM representation 29|]. Particular attention 
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has been devoted to the spectral properties and to the concept of SDs, which is expected 
to play an important role in explicit calculations also in nonextensive quantum many-body 
problems. Besides, we have presented the q EMM and the nonextensive version of the less 
known SDM j?] for a direct calculation of the q GFs and q SDs, respectively. A remark- 
able feature is that these methods should allow one, at least in principle, to explore on the 
same footing the nonextensivity effects for a wide variety of quantum many-body systems 
overcoming the a priori knowledge of the q partition function and hence of the q free en- 
ergy. Unfortunately, in contrast to the extensive case , the q CPs of two generic 
operators A and B cannot be expressed in a simple way, as the q GFs, in terms of the q 
SD KqAsi'^) defined in terms of the same operators. This introduces in the formalism some 
intrinsic difficulties which, in view of the present still limited experience, can be overcome 
only under appropriate constraints such as, for instance, close to the extensive regime. Un- 
der the restrictive condition (156|) . the calculation of q CPs is sensibly simplified and we have 
shown how the formalism works in exploring the g-induced nonextensivity effects on the 
low-temperature properties of the model (78)-(80) for the high-density Bose gas with strong 
attraction between the particles 31|. 

In any case, it is desirable to test again the effectiveness of the q EM and q SD methods for 
other many-body problems in nonextensive quantum statistical mechanics. A preliminary 



461] where we have applied the q SDM to an 



study along this direction can be found in Ref. 
isotropic Heisenberg model with long-range exchange interactions. 

We wish to stress again that, in the context of practical calculations, the crucial prob- 
lem remains to understand how one can calculate the relevant gCFs and the related q- 
thermodynamic quantities by using the general formulas (15^ and (!55l) beyond the limiting 
condition fl56l) . 

In conclusion, in light of t he g reat experience acquired in the extensive statistical me- 
chanics [ll, 2I, y, y, y, J 3, y, y, 10| , we beHeve that the q GFs technique and, in particular, 
the q SDM may constitute a powerful tool of investigation also in nonextensive many-body 
theory. 
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APPENDIX A: TWO-TIME GREEN'S FUNCTIONS AND SPECTRAL DEN- 
SITY METHOD IN NONEXTENSIVE CLASSICAL THERMOSTATISTICS 



As mentioned in the Introduction, the pioneering framework of the two-time GF method 
in extensive classical statistical mechanics by Bogoliubov and Sadovnikov Sj has opened 
the concrete possibility to describe classical and quantum many-body systems on the same 
footing. Besides, in many physical situations (when the quantum effects are negligible), 
the use of the classical formalism may offer substantial advantages especially from the com- 
putational point of view because in the calculations one handles only functions and not 
operators. 

Recently, the two-time GF technique and SDM have been formulated in nonextensive 
classical thermostatistics, within the OLM framework, in two our papers 



27 



and conve- 



niently applied to the ferromagnetic Heisenberg chain. For completeness, in this appendix, 
we review the basic equations to underline the main differences with the corresponding 
quantum case which has been the subject of the present article. 

First, we note that the basic ingredients of the quantum Tsallis thermostatistics shortly 
reviewed in Sec. [TTl remain formally valid for the classical framework. Here, p = p(q, p) 
denotes the probability distribution defined in the space phase F, q = {qi, qj\f} and p = 
{pi, ...,pj^} are the generalized coordinates, and Tr(...) stands for J {...)dT = J Hi^i dqidpi, 
for a classical system with Hamiltonian i/(q, p) and A/" degrees of freedom. As for the 
quantum case, also in this appendix we adopt the classical OLM representation [see Eqs. 

(HOD- (HI])]. 



1. Classical two-time q Green's functions and g-spectral density 

In classical nonextensive thermostatistics, the two-time q GFs involving two classical 



observables are defined by 
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G%{t,t') = {{A{ty,B{t')))['^ =eAt-t'){{A{t),B{t')})^, (Al) 

where the set of Poisson brackets {...,...} replaces the quantum commutator — z[. ...... ]^ 

and the q expectation value is taken over the phase space within the OLM spirit. Here, 
X(t) (with X = A,B) represents a dynamical variable depending on time through the 
canonical coordinates with X{t) = X(q(t),p(t)) = X(e*-^«*q(0), e*-^«*p(0)) = e*-^«*X(0), 
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where Lh = i {H, ...} is the Liouville operator. Hence, e*^^* acts as a classical time-evolution 
operator which transforms the dynamical variable X{0) = X at initial time t = into the 
variable X{t) at arbitrary time t, satisfying the Liouville equation dX{t)/dt = iLuXit). 

The main difference of the quantum and classical q GF methods lies in the relation 
between the q GFs and the q CPs. Actually, starting from the definitions of Poisson's 
brackets and q mean value, it can be shown that \^!^ 



GfUt, t') = qpe,{t - t')^F,AMt^ t' 



, d 



dt 



connecting G^^J^Q{t,t') to the new generalized two-time g-CP 



qAqB, 



{t,t') = {A,{t)B,{t'))^ 



(A2) 



(A3) 



of the two g-dynamical variables Aq and Bq defined (working in the canonical ensemble) by: 

X 



Cq = l-l3{l-q){H-Uq)>0. 



(A4) 



Unfortunately, it is not possible to derive a direct relation between G^^^^it, t') and FgAsit, t'] 
as happens in the extensive classical case ij]. Nevertheless, using the series representation 



^ CO 

= 1 + 2Vj„ 

1 — X ^ 

n=l 



(A5) 



where Jr, (z) are the Bessel functions 



\ 2m+n 



m=0 



m!P (n + m + 11 



(A6) 



with X = {1 — q)(3{H — Uq), we can formally write the expansion 



n,m=l 



where 



C7'"(/?) 



2(-l)™[f/3(l-g)]2'^+'^ 



(A7) 



(A8) 



m\V{n + m+ 1) 

With the restriction |x| < 1 the simplest series (1 — x)~^ = X^^o^" '^^^ used to obtain 
the formula 



G%{t,t')=q[3eAt-t') 



d 

dt 



{Amit'))q + $^r(l - qr {A{t)B{t'){H - UqT)^ 



n=l 



■ (A9) 
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In particular, we have 



d 



G^IbM ^ QPO^it - t')j^{A{t)B{t')),, (AlO) 



under the condition |x| = — q){H — Uq)\ <^ 1. This shows that only in the lowest order 
may the q GFs be directly related to FqAsit, t') as in the extensive case jS]. Any way, we can 
improve systematically the approximation ( lAlOl) . taking into account the successive terms 
in the expansion ( lAOp relating q GF to q CPs of increasing order. 

Assuming time-translational invariance, all the spectral properties derived in Sec. IIII Al 
remain formally unchanged, but now the classical q SDAg^^ (u) is defined by 

/ + 00 
dTe'-^{{A{T),B})^. (All) 
-oo 



Besides, from Eqs. (IA2p and flA3p . one can show that |27| 



(^WB,>.= /'°°^^^^2M£_. (A12) 

oo 



2n qPu 

Then, using the series representation ( lA9p . we get 



J-oo qP n,m=l 



which corresponds to Eq. ( IMl) and allows us to connect the classical q GFs with the related 
q SD. In the limit g — > 1, Eq. flAlSp reproduces consistently the well-known exact extensive 
result [9 1 

lu„(AB), = (AB)= rf^^H£M. (A14) 

In order to derive the spectral decomposition for KqAB (^^), we now introduce a Hilbert space 



Sq of the classical dynamical variables with a scalar product defined conveniently by 
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{A\B)q = Zq{A*B)q, (A15) 

with Zq = J dT[l - I3{1 - q){H{q,p) - Uq)]—^. In this space one can consider the eigen- 
value equation Lh'^^ = ^k^k for the Hermitian Liouville operator. It is also immediate 
to prove that, if is an eigenfunction of with (real) eigenvalues tUfc, then is also 
an eigenfunction of with eigenvalue —uJk- If we assume that {^k} is a complete set of 
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orthonormal eigenfunctions, we can consider the expansions A{q,p) = yii.(^fc|^)g^fc(q, p) 
and _B(q, p) = ^^(\E'fc|i?)q\I/fc(q, p). Bearing this in mind, we can write [27 1 



A,ab{oo) = 2'Kq(3ujZ^^ Y,{mk\B,),{^l\A,),6{uj - uu), (A16) 

k 

which is the desired classical g-spectral decomposition for KgABil^)- particular, ii B = A* , 
we have that KqAA*i!^) is a real and positive-definite quantity. A consequence of Eq. (1A16P 
is that GqAsi!^) {Ag{t)Bg)g can be written as 



GgABiu^) = 2nqPZ^' ^(^,|s^)^(vl>* (A17) 

to — UJk 



^[^k\-t^q)q{'i^k\^q)q-^'' 

k 

and 

{Ag{T)Bg)g = Z " ^ ( , | 5, ) , ( V[/ ^ | A, ) , 6 ^ ^ ^ (AlS) 

k 

Thus, also in the classical case, the real poles of GqAsi'^) represent the frequency spectrum 
of undamped oscillations. 



2. Methods of calculation for classical two-time q Green's functions and g-spectral 
density 

a. Classical equations- of -motion method 

As in the quantum case, successive differentiations of Eq. (lAip with respect to r = t — t', 
yield the infinite hierarchies of classical coupled EMs in r and uj representations: 

^ {{C^A{t)- B))'^^^ = 6 (r) B})^ + {{C-^'A{t)- B))^^^ (m = 0, 1, 2, ...) (A19) 

and 

u: {{C^Air); B))^^l = i B})^ + i {{CT'A{r)- B)f^l (m = 0, 1, 2, ...), (A20) 

respectively, where Ch = iLn and C]}A means C^^A = A, C\jA = {A, H}, Cj^A = 
{{A, H} , H} and so on. 

At this stage, the considerations made in Sec. IIVBI for solving the EMs in the quantum 
case apply for the classical one, too. 
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b. Classical q-spectral density method 



From the definition Ag^B (t) = i {{A{t), B})^ of the classical q SD in the r space, suc- 
cessive derivatives with respect to r yield 



AgABir)=i{{C'SA{T),B})^ (m = 0,1,2,...). (A21) 



Then, proceeding as in the quantum case, we easily find 

/ + 00 T 
— a;"^A,^s M = -{^r-' {{C'SA, B})^ (m = 0, 1, 2, ...). (A22) 

The quantity on the left-hand side will be called the m moment of Ag^B {^^), and relation 
(1A22I) can be seen as an infinite set of exact MEs or sum rules for the classical q SD. 

As for the quantum counterpart, due to the possibility to evaluate the Poisson brackets 
and hence the q averages on the right-hand side of Eq. flA22p . the m moments can be explic- 
itly calculated without a priori knowledge of AgAB (t^)- So one can consider the sequence 
of equations, Eq. flA22p . as a typical moment problem to determine the unknown function 



AgAB (i^) and hence all the related macroscopic ^quantities. At this stage, in view of the 
classical exact g-spectral decomposition ( ]A16[) for AgAB {^), the basic idea of the SDM [7] for 
classical gthermostatistics does not differ from the quantum case and all the considerations 
made in Sec. IIVI about the possible functional representations for the q SD preserve their 
validity in the context of the classical q many-body theory. 

APPENDIX B: CALCULATION OF k SUMS AS F ^ oo IN THE HIGH- 
DENSITY AND LOW-TEMPERATURE LIMITS 

For the (T = 0)-sums over k in Eqs. ffTTTD-ffTT3D. with l^"^ Ek^o(--) = 
{2tt'^)^^ dkk'^{...) as V ^ oo (all the summands depend only on A; = |k|) and the trans- 
formation X = k{n'j)~^^^, we can write 

1 H-fi^f 



^yj}L^ = y!n2i dx ^' ; ^ , (bi) 



(B2) 
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2V u;w (27r)' 



where 



fix) 



00 X 

dx — 



1 2 







(n7) 2 



-(n7)i/2^2fl! 



(B3) 



(B4) 



(727)^2 + ?7 + 1 

Following Babichenko 3l|, with {wyY^"^ ^ 1 one can now neglect the term (717)"^/^ with 
respect to and the Gaussian part in Eq. flB4l) . This is quite legitimate since the basic 
contribution to the previous integrals is made by x ~ 1. 
Then, one finds 

/k - t^k 



k^O 



a(n'y) 



^5^(/k-^k) ^ ^^(^7)^ 



kT^O 



where 



1 



(27r) JO 
1 



kT^O 



dx 



5 



(B5) 
(B6) 
(B7) 

(B8) 

(27rrJo + 477^/' ^4' ^^^^ 

Concerning the T-dependent sums over k in Eqs. f lllll) - flll3p . in the low-temperature limit 
only small values of k contribute to the integrals which contain exponential functions. So, 
with the additional condition (717)^/^ ^ 1, all the integrals can be exactly computed and we 
get, to leading order in T, 



487rV2 ^4^ 



-E 



V ^ _ 1 

kT^O 



T/ 



/k 



30 V T] 
1 + l 



12 



E 



27r2 /r7 + 1 



k^O 



V 



E 

k^O 



1)' 



15 V ?7 
1 /r/ + 1 



(e/3^k _ 1)^ 



6 



1 
3 

'(n7)-tT^ 



(BIO) 

(Bll) 
(B12) 
(B13) 
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1 ^ huje^^^ _ 27r2 

Inserting the previous results into Eqs. ( ]llip -( fTT3l) . one immediately obtains expressions 
flll4l) - flll6p for the depletion of the condensate, the chemical potential, the q internal energy 
density, and hence the q specific heat flll7p as a function of T and n in the high-density and 
low-temperature limits. 



(B14) 



r/ + 1 



(B15) 



{n-i)-^T^. 



(B16) 
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